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The dynamical mean-field theory together with the non- 
crossing approximation is used to set up a novel scheme to 
study the electronic structure of strongly correlated electron 
systems. The non-interacting band structure is obtained from 
a density functional calculation within the local density ap- 
proximation. With this method the doped Mott insulator 
Lai-xSrxTiOs is studied. Starting from first-principle calcu- 
lations for a cubic and an orthorhombic system we determined 
the one-particle spectrum. Both one-particle spectra show a 
lower Hubbard band (seen as d° transitions in photo 

emission experiments), a Kondo resonance near the Fermi en- 
ergy and the upper Hubbard band (d^ transitions in an 
inverse photoemission experiment). The upper Hubbard band 
develops a multipeak structure, a consequence of the consid- 
eration of all local two-particle correlations, which leads to 
the full multiplet structure in the atomic limit. The calcu- 
lation for the orthorhombic system shows qualitative good 
agreement when compared with experimental photoemission 
spectra. 

71.27.-|-a Strongly correlated electron systems , 71.30.-|-h 
Metal-insulator transitions and other electronic transitions , 
74.25.Jb Electronic structure 

I. INTRODUCTION 

One of the most challenging problems in condensed- 
matter physics is the accurate calculation of the elec- 
tronic structure starting from first principles, since it is 
impossible to solve the many-body problem without se- 
vere approximations. A great success in first principle 
methods was ihe development of the density functional 
theory (DFT)EJ, which produces an exact groundstate en- 
ergy and maps the many-body problem onto a ficticious 
single-electron problem with an one-electron exchange 
and correlation potential. The method would be exact, 
if one knew this particular potential. 

A commonly used approximation replaces this func- 
tional by the one of the homogeneous electron gas. This 
approximation is widely knowri as the local (spin-) den- 
sity approximation (L(S)DA)El. The LDA turns out to 
be very crude, if the on-site Coulomb interaction be- 
tween electrons becomes very strong, i.e. if the interac- 
tion strength is of the order of, or larger than the band- 
width. This is particular the case in systems with d- or 
/-electrons. As is well known, the L(S)DA approxima- 
tion leads to wrong results for stronlgy correlated elec- 
tron systems such as Mott-Hubbard insulators or heavy 



fermion systems. 

A first step towards a better description of strongly 
correlated systems was the introduction of the LDA-I-U 
methodcl, a Hartee-Fock-like scheme, that is able to ac- 
count for a variety of interesting effects in transition 
metal compounds. For example, the LDA-I-U allows 
to describe correctly the groundstates of magnetically, 
charge and orbitally ordered materials. Due to its 
Hartree-Fock nature, however, this approximation still 
has severe drawbacks. 

To gO|-lDeyond LDA-I-U, a new scheme was recently 
proposeda, which combines the LDA bandstructure cal- 
culation with the dynamical mean-field theory (DMFT: 
for a review see reference 5). The DMFT maps the lat- 
tice model onto an effective impurity problem and can be 
used to solve the many-body problem self-consistently 
with various methods, like e.g. quantum Monte Carlo 
(QMC)eI, exact diaspnalization (ED)td, iterated pertur- 
bation theory (I|PT)l3 or within the non-crossing approx- 
imation (NCA)EI. 

A calculation scheme for LDA+DMFT in conjection 
with IPT was already accomplished in reference 4. How- 
ever, in this work only one mean Coulomb parameter U 
was used to describe all local correlations. It is however 
well known that the complexity of the physics observed 
in e.g. transition metal oxides is governed by the subtle 
interplay of the various local Coulomb integrals, i.e. the 
neglect of the detailed orbital structure presents a serious 
limitation of the method. 

In this paper a LDA-I-DMFT approach with the NCA 
is introduced, which allows to consider the full local mul- 
tiplet structure. Thus, an intra orbital Coulomb energy 
U, an inter-orbital Coulomb energy U', and Hund's cou- 
pling constant J are considered within a SU(N) invari- 
ant formulation of the local Hamiltonian for degener- 
ated bandsEJ. The bandstructure calculations as starting 
point for the DMFT were dcpe with the linearized muffin- 
tin orbital method (LMT0)E2l. The paper is organised as 
follows. In sectional we introduce the Hamiltonian used 
for the many-body calculation and describe the details of 
the NCA scheme used to solve the local problem . Results 
for Lai-xSrxTiOa as a particularly interesting example 
of a correlated me tal obtained with this technique are 
collected in section HI and a summary in section ^ will 
conclude the paper. 
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II. THE CALCULATING SCHEME 

In order to combine the achievements of the DMFT 
with the framework of DFT, one needs a calculation 
scheme that allows to map the calculated band struc- 
ture onto a suitable tight-binding model. This can be 
done most conveniently with the LMTO in the orthogo- 
nal approximationllil. With this technique a tight-binding 
model can naturally be constructed from the DFT+LDA 
bands in real-space representation 
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Since the LDA one-electron potential is orbital indepen- 
dent and takes into account the Coulomb interaction in 
an averaged way, one can easily generalize this Hamilto- 
nian by adding local Coulomb correlations: 
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In term (2) and (3), the index i labels the unit cell, / de- 
notes the atom in the unit cell and ma and others stand 
for orbital and spin. The prime on the sum indicates that 
at least two of the indices on different operators have to 
be different to account for Pauli's principle. In our ap- 
proach we are able to include all local two-parcticle corre- 
lations appearing in term (2) and (3). In the following we 
will assume, that it is only necessary to take into account 
the Coulomb interactions for the d-shell of the transition 
metal ions {i = id and / = Id) explicitly. Therefore the 
indices il will be omitted. All other orbitals will be con- 
sidered as itinerant bands, which are well described by 
the LDA. 

Let us discuss the interaction term in more detail. The 
first term (2) describes the density-density Coulomb rep- 
sulsion. Here, we introduce two distinct Coulomb pa- 
rameters: the intra-orbital Coulomb energy U has to be 
considered in case of a doubly occupied orbital, while the 
inter-orbital Coulomb energy U' applies for example in 
the case of a doubly occupied d-shell with electrons on 
d-orbitals with different indices. The second term (3) de- 
scribes Hund's coupling and all other off-diagonal inter- 
actions. Although in principle the full Coulomb matrix 
(2) and (3) can be used, we restrict ourselves for practical 
purposes to the following, commonly used form for the 
two-particle interactions of the d-electrons 
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with a density-like Coulomb term (4), a SU(2)-invariant 
form of Hund's Coupling (5) and a charge-flip term or 
pair-hopping term (6). Jc — J '^s chosen here as a good 
approximation. 

Since the LDA already contains the influence of the 
Coulomb interaction to a certain degree, the problem of 
double counting of these contributions by using Hcorr 
arises. In order to avoid this double counting one has to 
subtract off the interaction contributions from the LDA. 
Unfortunately, the precise form for one particular set of 
orbitals is not known and the best one can do is to ac- 
count for these contributions in an averaged way. Since 
the LDA total energy is to a good approximation a func- 
tion of the total number of electrons, one can assume that 
the interaction contributions have the form 



Ei = ^Und{nd 
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Here, U is the mean value for the Coulomb interactioii 
and may be obtained from a first-principle calculationEHI 
or from experiment, for example high-energy spec- 
troscopy. Ud is the total number of d-electrons and nd,a 
is the number of c?-electrons with spin a. Distinguishing 
intra- and inter-orbital interaction for a iVdeg-fold degen- 
erate orbital-system one can determine U and U' by the 
relation 



U = 



U + 2{Ndeg - l)U' 
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(8) 



Starting from given values for U and J and using the rela- 
tion U = U' + 2J, which is a coiKequence of the rotational 
invariance of the HamiltonianliJ, one is able to determine 
the intra- and inter-orbital Coulomb interaction U and 
U'. Finally, only an expression for the local energy of 
the d-states is needed. In the spirit of L(S)DA-t-U we 
can define such a one-electron energy by 
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and Elda the total energy as calculated from 
DFT+LDA. With the corrected total energy functional 
Elda — E] one has to construct the tight-binding Kaniil- 
tonian within the framework of the LMTO-methodlld. 

Given the band structure we are now in the position 
to set up the scheme necessary for the DMFT. First, let 
us define the non-interacting Green function 



G°(z) = / diu 



(12) 



via the Hilbert transform of the spectral function 
obtained from the Hamiltonian (1). Note that in gen- 
eral PctI"-^) ^^'^ consequently Gj^(z) will be matrices in 
orbital space. The most important feature of the DMFT 
is that the proper one-particle self esergy due to the local 
Coulomb interaction is purely locaB. Thus, we obtain as 
an expression for the full Green function of the interact- 
ing system 



G,(z) = G°(z-E(z))= / duj 



(13) 



This equation must again be read as a matrix relation in 
orbital space. In the spirit of the DMFTQ eq. (13) can 
be cast into the form 
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where a so-called hybridization function was introduced, 
which fulfills the standard relation 
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The one-electron energy e^, the imaginary part of the 
hybridization function and the local Coulomb parameters 
define a multi-orbital impurity model, which has to be 
solved to obtain the selfenergy ^a{z). Finally, eqs. (13), 
(14) and the solution of this impurity model define a 
self-consistency cycle for the solution of the lattice model 
within the DMFT. 

In the following we will set up a resolvent perturba- 
tion theory for the multi-orbital impurity model defined 
above. Within this approach the model is solved approx- 
imately by considering the lowest order diagrams only, 
which is the well-known NCAlI. To this end we first have 
to diagonalize the local Hamiltonian 



Hioc = ^Ea\a){a\ 
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and use the basis states {\a)} in the further procedure. 
As a consequence, all local interactions are fully included 
in these impurity states. Therefore, all types of local 
interaction terms can in principle be considered leading 
to a realistic multiplet structure in the atomic limit. To 
account for the hybridization term, we have to express 
the fermionic creation and annihilation operators for the 
d-orbitals in terms of our basis {\a >}, i.e. 



Q,/3 



(17) 
(18) 



with K the orbital index and cr for the spin. With (17) 
and (18) |«ne can straightforwardly formulate the NCA 
equationaa as 



i?„(z) = i?^(z) + i?° (z)E„(z)i?„(z), 



with 
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/(e) is the Fermi function and rjaa' is equal to 1 or - 
1 depending on wether one has to add or subtract an 
electron to excite the impurity from state la) to state \a'). 
Again, eqs. (19)-(21) have to be solved self-consistently. 
Finally, we are able calculate the c?-Green function from 
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In eq. (22), Z = 



dze->^' 
2111 



Ra{z) denotes the impu- 



rity partition function and (3 is the inverse temperature. 

The solution eq. (22) of the impurity problem closes 
the cycle (13) and (14): Starting from an initial guess 
for the selfenergy T,cr{z) one obtains from eq. (14) the 
hybridization function A^(z), which is inserted into the 
NCA (21). From the impurity Green function (22) and 
the hybridization function one can then calculate a new 
self energy from eq. (14) and repeats the procedure until 
self consistency is reached. 



III. RESULTS FOR LAi-xSRxTIOs 

We applied the procedure discussed in the previous 
section to Lai-xSrxTiOs. This material is at high tem- 
peratures a paramagnetic metal for the undoped case 
(x=0) and becomes an antiferromagnetic insulator below 
Tn — 125K with a very small gap of 0.2eV. Doping with 
a few percent of Sr leads to a paramagnetic metal with 
a large effective masalj. As photoemission experiments 
show, Lai_xSrxTi03 may be regarded as an example of 
a strongly correlated metal, i.e. the density of states is 
that of a doped Mott-Hubbard insulator. 
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The LaTiOa crystal is a slightly distorted cubic per- 
ovskite. The octahedral coordination of the oxygen ions 
leads to a t2g-eg crystal-field splitting, such that both 
bands t2g and eg are well separated, as can be seen in 
a standard LDA calculation (LMTO method). At the 
Fermi level one finds the t2g band of the Ti id states. 
3eV above the Fermi level and well separated from the 
t2g band the Cg band is located. Lead by these results we 
can restrict ourselves to the 3-fold degenerate t2g orbitals 
as relevant for the low temperature properties. 

The Coulomb parameters U — AeV, J — O.GeV we use 
for the DMFT calculation were obtained from a con- 
straint LDA calculatior£3. U was determined in a calcu- 
lation, where all electrons except t2g ones can screen the 
Coulomb interaction. ( The influence of the Cg electrons 
on the CcuilaEpb interaction inside t2g subshell was inves- 
tigated inEj'Eil). In order to show the effect of varying the 
Coulomb parameter on the results we chose U = 6eV in 
the first parameter sets, which refers to the case where 
Cg-electrons do not participate in the screening. 



cubic 




S) in [eV] 

FIG. 1. The non-interacting t2g 

density of states for LaTiOs corrected according to eq. (7) 
with U = 6eV, J — O.QeV. The upper part shows the density 
for a cubic, the lower part for an orthorhombic system. 

In figure |l| we show the non-interacting density of 
states for the t2g bands only taken from a LDA calcula- 
tion for a cubic (upper part of figure |l|) and a orthorhom- 
bic system (lower part of |^). Both spectra are corrected 
according to eq. (7). From the DOS in figure |l| we can 
determine the one-particle energy e'}^^ as average energy, 
and the hybridization function At2g{uj)is then defined 
by equation (13) and (14). In figure ^ the bandwidth 
determined by the main contributions of the spectrum is 
varying for different underlying structures examined by 
the LDA calculation. The orthorhombic system leads to 
a narrower spectrum with the main contribution below 
the Fermi energy. The spectrum for the cubic system is 
broader and the spectral weight is spread over the whole 



range of the band. As a consequence, the average energy 
of the orthorhombic system is lower than for the cubic, 
which leads to different one-particle energies in the re- 
sulting impurity model. It is evident that this change 
of bandwidth and of one-particle energy will eventually 
lead to different physics for both systems. 

A particularly important simplification arises from the 
fact that the tight-binding Hamiltonian (1) constructed 
from the corrected LDA t2g bands turns out to be diago- 
nal in orbital space. This means that all orbital matrices 
in eqs. (13), (14), (21) and (22) are diagonal and - since 
we do not allow for orbital ordering ~ degenerate with re- 
spect to the orbital indices, which results in a tremendous 
reduction of the necessary computational effort. 
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FIG. 2. t2g density of states for Lai_2,Sra;Ti03 for cu- 
bic and orthorhombic crystal structure. The parameters 
are: U — 6eV and J = O.GeV, a reciprocal temperature 
f3 — 100 eV~^ and a doping of x = 6%. The denoted tran- 
sitions are the removal of electron from a single occupied 
t2g-state and the addition of an electron to a doubly occu- 
pied state. 

The results of the self-consistent DMFT calculation for 
the t2g spectral function of a cubic and a orthorombic sys- 
tem with a partially screened mean Coulomb interaction 
of U = 6eV and J = 0.6eV and a reciprocal temperature 
(3 — 100 eV^^^ are shown in figure^. The doping level was 
set to x = 6%. Both systems lead to qualitatively similar 
spectra. Especially the lower Hubbard band reaches from 
low energies up to the Fermi level. This band refers to 
d'^ transitions in a photoemission experiment. At 
the Fermi energy the well known many-body resonance 
can be seen, which describes singlet-triplet excitatieirs of 
the impurity and the surrounding effective mediumEj. At 
higher energies the upper Hubbard band occurs which 
refers to d^ transitions in an inverse photoemis- 

sion experiment. Here contributions, which stem from 
transitions into several multiplet states, are almost com- 
pletely smeared out to a broad, featureless band. Apart 
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from these gross features, the spectra show clear differ- 
ences between orthorombic and cubic system. Especiahy 
the influence of different bandwidth is clearly visible, the 
narrower non-interacting spectrum for the orthorhombic 
system leading to narrower peaks in the DMFT spec- 
trum. As a consequence, the value of the gap between 
lower and upper Hubbard band is increased, too. Let 
us stress that the value of the gap, due to the complex 
atomic multiplet structure, is not simply given by the 
pure Coulomb interaction U. For a partially screened 
mean Coulomb value oi U = 6eV the real Coulomb pa- 
rameters are U = 6.96ey and U' — 5.76eV. Thus the 
lowest lying excitation is into a spin S = 1 state with 
an excitation energy of e°^^ — {U' — J). This would lead 
to a gap value of 5.l6eV minus bandwidth. The smaller 
bandwidth of the orthorhombic sytem thus turns out to 
produce an overall larger gap. 

In order to exploit the difference of the value, and 
hence the impact of the procedure to determine U, we 
compare the previous results to a calculation where a 
fully screened mean Coulomb energy of C7 = ieV is used 
(i.e. U = A.96eV and U' = 3.76eV). Note that the value 
of J is not changed by using different screening scenarios, 
which reflects the well-known fact that the exchange cou- 
pling is rather insensitive to the detailed orbital structure 
and chemical surrounding in transition metal oxides. 



interaction compared to the partially screened case. The 
Kondo resonance is enlarged as a fact of different one- 
electron energy, which is the lowest for the orthorhombic 
set and for the full screened Coulomb energy U. 
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FIG. 3. The t2g density of states for LaTiOs for the or- 
thorhombic structure and U = 4eV and J — 0.6eV at a re- 
ciprocal temperature /3 = 100 eV~^ . 

The rcsuh of the LDA-fDMFT calculation for the or- 
thorhombic system is shown in figure ^ The lowest en- 
ergy difference between and cP is now reduced by 2eV 
leading to a total shift of the upper Hubbard band com- 
pared with figure ||. The multiplet structure is not af- 
fected since it is proportional to the coupling constant J 
only. For this set of parameters at a doping of a; = 6% the 
lower and the upper Hubbard band is shifted towards the 
Fermi energy as a consequence of the smaller Coulomb 
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FIG. 4. The t2g density of states for LaTiOs for orthorhom- 
bic structure given by the LDA-I-DMFT calculation (full line) 
compared with LDA results (dotted curve) and data points 
from a photoemission epsperiment for Lai-xSrxTiOa at a dop- 
ing of a; = 6% (circles )Ej. The parameters of the calculation 
were U = 4eV and J = 0.6eV and a reciprocal temperature 
f3 = 100 eV-^ 

A comparison with the experiment is done in figure 
^ using the spectrum of the orthorombic system. The 
bandwidth, maximum and center of mass of the lower 
Hubbard band are in fair agreement with experiment. 
The many-body resonance at the Fermi level in the spec- 
trum in figure |2| manifests itself as an additional structure 
at the Fermi level as a consequence of the convolution 
with the Fermi function. Note that such a structure - al- 
though less pronounced - can bee seen in the experimen- 
tal data, too. There are, of course, still differences in the 
exact distribution of the spectral weight between exper- 
iment and our theory, which possibly originate from our 
approximation of a purely local self-energy. Nevertheless, 
compared to the pure LDA result, which gives a com- 
pletely wrong account of the spectrum, the LDA+DMFT 
appears to capture the essentials of the physics in this 
material. 
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IV. SUMMARY 

In this paper we have described one possible reaUza- 
tion of a combination of density-functional theory with 
the LDA approximation and the recently developed dy- 
namical mean field theory to obtain a first-principles 
calculation scheme for strongly correlated electron sys- 
tems. To solve the DMFT equations we used the Non- 
Crossing approximation, which allowed us to include all 
local Coulomb interactions in our calculations, which 
is important to reproduce the full multiplet structure 
in the atomic limit. In order to test the method, we 
showed results for the doped Mott-Hubbard insulator 
Lai_xSrxTi03. The comparison of our results with a 
photocmission experiment shows, that the lower Hub- 
bard band is of the correct width and the center of mass 
of both experiment and theory is at the same position. 
Also, a shoulder seen at the Fermi energy in experiment 
is accounted for by our results as the onset of a many- 
body resonance. In addition, the method appears to be 
very sensitive to the atomic parameters and is clearly 
able to discriminate between different crystal structures 
and symmetries, although the lattice enters the DMFT 
only in an averaged manner. In total the results are very 
encouraging for the perspective of future applications. 
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